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I Abstract. Using optical model calculations it has recently been shown that refractive phenomena 

. from the collective mean field can significantly alter the sizes inferred from two-pion correlations. 

^ ■ We demonstrate that such effects can be accounted for in classical calculations if mean field effects 

I are included. 
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^ ! The six-dimensional two-pion correlation function C(P, Q), measured as a function 

Q I of the total and relative momenta, has provided crucial insight into the space-time 

O ■ development of heavy-ion reactions [1]. Recently, the refractive effects of particles 

^ . traversing the mean field were calculated by symmetrizing quantum waves in the frame 

ly-^ I work a time-independent complex optical potential [2, 3]. If the pions left a region 

O ■ of lower mass, it was found that the lensing effect of the potential was to distort 

the extracted source dimensions by a few tens of percent. The effect was analagous 
T ! to lensing effects from Coulomb mean fields, which were shown to explain different 

Q I apparent source sizes for positive and negative pions at AGS energies [4, 5]. 

^ ■ The refractive corrections studied in Reference [2, 3] allowed for a more physical 

', interpretation of correlation data from RHIC. When ignoring the corrections it appears 

^ I that the fireball grows to a transverse (to the beam) radius of ~ 13 fm, expanding at a 

^ ■ speed ^ 0.7c, and rapidly disintegrates at a time of ~ 10 fm/c [6]. The puzzling aspect 

of this picture is that the fireball surface must expand at a speed of 0.7c from the initial 
collision to grow from its initial size of 6 fm to 13 fm in 10 fm/c, not allowing any time 
for the matter to accelerate. Correlation analyses typically provide three dimensions: 
R\ong, the longitudinal size defined parallel to the beam, i?out, the dimension of the phase 
space packet perpendicular to the beam and outward along the direction of the pair's 
momentum, and R^ide, the sideward dimension which is perpendicular to both the pair 
momentum and the beam axis. The apparent sideward dimension in [2, 3] was shown to 
be increased by the refractive effects of the potential, which would suggest the true radius 
of the fireball was somewhat smaller than the 13 fm previously believed. Reducing the 
size by one or two fm would provide a more physically plausible picture of the reaction's 
evolution. 

The calculations in Ref . [2, 3] involved solving for single-particle outgoing wave func- 
tions in the presence of complex optical potentials, then using the symmetrized product 
as a correlation weight. Quantum calculations such as these have some drawbacks. First, 
implementation becomes complictated if one were to account for time- varying poten- 



tials, or potentials without the spherical and boost symmetries assumed in [2, 3]. Sec- 
ondly, finding expressions for the mean field, especially the imaginary part of the optical 
potential, would be quite involved. Finally, it is difficult to extend such calculations be- 
yond the case of identical particles as one then needs to consider the quantum three-body 
problem. For these reasons, we wish to study whether the effects can be calculated by 
considering classical trajectories through the mean field, then using the asymptotic phase 
space density to generate correlations. Given that the inferred diameters of the RHIC 
fireballs are often near 25 fm, one might expect that classical considerations could be 
valid except at very low pt. 

For non-interacting identical particles, two-particle probabilities can be calculated in 
terms of the outgoing phase space density [1]. 



C(P,Q) = l + y'jV^p(r)cos(Q'-r), (1) 
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where the primes denote the positions measured in the pair frame, and Q' = — v't 
is the relative momentum in that frame. Since describes the separations of particles 
with the same velocity, the expression is independent of t' once it satistifies the criteria 
of being beyond the point at which particles interact with the remainder of the system, 
including mean-field interactions. 

An alternative approach is to treat interaction with the mean field separately, define the 
source function so that it describes the points at which particles had their last interaction 
other than through the mean field, then use the outgoing wave function, which is a 
solution to the equations of motion using the mean field, to describe the evolution from 
X to the asymptotic momentum state [2, 3, 4]. This involves replacing the phase factor 
used to describe the evolution from x to its asymptotic state p. 

The two-particle probability is then, 
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The correlation function is the ratio of the two terms. Here, is the outgoing single- 
particle wave function describing evolution through the mean field. The source functions 
s{p,x) now refer to the points where a particle had its last non-mean-field interaction. 

Both approaches should be equivalent if the classical trajectories would correctly de- 
scribe propagation through the mean field. We present an "apples-to-apples" comparison 
of the two approaches using the same complex optical potential, 

U{r)^{U, + iUj)-^^-^^, (4) 



where a is the diffuseness parameter. If a = 0, the potential is a step function. We solve 
for the distortion to the apparent sideward radii for particles which originate isotropically 
from points a radius ro from the center of the cylinder. For identical particles, radii can 
be determined from correlations using the expression [1], 
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where the radii represent the dimensions of the asymptotic phase space density, not the 
dimensions of the emission points. Using Eq. (3) which gives the correlation function 
in terms of the distributions points for last collisions convoluted with outgoing wave 
functions, one can then apply Eq. (5) to write an expression for the variance of the 
sideward size, 

where the asymptotic momentum p moves along the x axis and are asymptotic 
outgoing wave functions. 

Calculations were performed for s being independent of the direction of p with 
emission points being confined to a radius of tq. The potential parameters were Ur — 
'"med ~ '"vac ^^^h the in-medium mass mmed = 50 MeV/c^ and Ui = ■ Tq with 
To = 100 MeV. This corresponds to a classical decay rate of F = Yqih-h/E. The fact 
that the decay rate falls ~ 1/i? is characteristic of a scalar form for the optical potential. 
The lower panel of Fig. 1 shows the distortion of inside- defined as the ratio of i?side with 
the potential to R^{^^ with ?7 = 0, as a function of pt. The distance scales for the potential 
were R— 10 fm and a = 3 fm. Results are shown for three different values of ro: 5 fm, 
10 fm and 15 fm. The distortions rise at low pt with stronger distortions for small tq. 
Also displayed in Fig. 1 are escape probabilities. Quantum mechanically, the probability 
that a particle escapes the fireball is 
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Since low pt particles spend more time in the fireball, and since the decay rate falls as 
l/E, escape rates are small at low pt. 

The corresponding classical calculations were performed by solving for the trajec- 
tories of particles through the mean field and are shown alongside the corresponding 
quantum calculations in Fig. 1. The apparent sideward source sizes from the calcula- 
tions used for Fig. 1 agree well with one another, differing by only a percent or two for 
/7j > 100 MeV/c. Even at pt = 40 MeV/c, the calculations agree to better than 10%. As 
explained above, the agreement should be expected to worsen for smaller a. Repeating 
the calculations for a = 0.5 fm, discrepancies increased to the level of a few percent for 
Pt > 100 MeV/c, and notably higher for pt ~ 50 MeV/c. However, it is hard to physically 
motivate such a sudden change in densities. 

One can also understand the connection between classical and eikonal approaches 
such as those discussed in Ref.s [7, 8, 9]. Since eikonal phases appear in interferences 
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FIGURE 1. The ratio of the apparent sideward dimension to the dimension without mean field is 
shown as a function of pi for the optical potential described in the text. Distortions are calculated for 
both quantum calculations (circles) and classical trajectory calculations (squares). The source function 
describing the final collision points was confined to an intial radius ro=5,10 or 15 fm. The upper panel 
shows the probability that such particles escape without being absorbed due to the imaginary part of the 
optical potential. The two approaches agree within a few percent. 



as d{pi,x) — 8{p2-,x), and since the derivatives of phase shifts can be identified as the 
spatial offset incurred by slowing down or speeding up in a classical path through the 
potential, an connection between eikonal and classical approaches ensues for pi ^ p2 
[10]. 

In order to better illustrate the physics of the refractive distortion we consider a simple 
example where pions are emitted from the surface of a cylinder of radius R, where the 
in-medium mass of the pion at the surface is m^^^ and the asymptotic vacuum mass is 
mvac- We assume the mass returns to its vacuum value exponentially, 

m\r) = ml, + (m^d - mlJe-^^-^y^. (8) 

From time-reversal arguments a thermal source will emit particles of given momentum 
with the same trajectories as those that describe absorption. The asymptotic trajectories 
of particles with momentum px = pt, Py = 0, that intersect with the cylinder have a 
uniform distribution of impact parameters up to a maximum b as illustrated in Fig. 2. 
The effect of an attractive mean field is to stretch the width of the phase space cloud 
by a factor bmax/R- The sideward dimension measured in two-particle correlations is 
stretched by this factor. 
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FIGURE 2. Trajectories that lead to the same asymptotic momentum are altered by an attractive mean 
field (left). The phase space distribution for particles of this given momentum are widened in coordinate 
space by a factor bmax/R- Repulsive mean fields (right) reduce the size of the region of the outgoing 
particles. 
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FIGURE 3. Assuming a static cylindrical source, the distortion to the sideward dimension due to the 
mean field is shown for attractive scalar fields. The field lowers the pion mass to 50 MeV/c^ and falls off 
exponentially outside the emitting radius R with a distance scale a. The distortion is stronger for larger 
ranges a and at lower pf . 



The maximum impact parameter for capture can be found by combining conservation 
of energy and angular momentum, L = bpt, 

pj = pj + V,ff{r) (9) 

^2 2 

One can analytically solve for the maximum impact parameter Z^max from which one can 
overcome the potential barrier with a give pt. The high pt behavior is independent of a, 
and in the case where a> R/2, the result is independent of a for all pt- 

Figure 3 shows the ratio bmax/R as a function of pr for several values of a given the 
case of a lighter in-medium pion mass, mmed = 50 MeV/c^. For a = 0, the mean field 
has no effect as the capture cross section does not extend beyond the cylinder. For the 
saturating value, a= R/2, the cross section is significantly enhanced. For a = R/IO, the 



enhancement is identical to the saturating value at large pt and is somewhat reduced at 
low Pt. The enhancement of the apparent sideward size is characteristic of an attractive 
mean field. This is consistent with Liouvilles theorem, which states that contraction of 
the phase space density in momentum space as it leaves the attractive mean field should 
be accompanied by a growth of the phase space cloud in coordinates space. 

The simple cylindrical picture was expanded to include both longitudinal and collec- 
tive flow [10]. For emission during the expansion stage, distortions were increased in 
magnitude since the boundary layer expands behind the pions, staying in the region of 
non-zero mean field for a longer time. Pions emitted during the collapse of the breakup 
surface were less affected since they spent less time near the surface. 

The calculations presented here were only meant to be illustrative. The potential 
magnitude of mean-field effects and the ability of classical pictures to model the effects 
underscore the importance of incorporating mean field effects into hadronic Boltzmann 
descriptions of the breakup stage. Such calculations have been undertaken in the past, 
especially for lower-energy Fermi-velocity collisions, where the mean field is paramount 
[11, 12]. In this energy range, phase space points used to generate correlations are 
typically taken from when a particle leaves the region of mean field, rather than when 
they have their last collision. Given the likelihood that the dispersion relations for 
pions are non-trivial, it is important that calculations incorporate momentum-dependent 
interactions. Such an effort is crucially important for reconstructing the space-time 
evolution of the fireball, and for understanding the pressure in the hadronic phase. 
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